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Abstract 

In  this  paper  we  present  an  analytical  evaluation  of  the  perturbations  to  mesoscale  atmo¬ 
spheric  flows  induced  by  thermal  inhomogeneities  in  the  convective  boundary  layer.  We  study 
the  time  evolution  of  these  perturbations  as  a  function  of  the  intensity  and  of  the  horizontal 
and  vertical  scales  of  the  diabatic  forcing.  The  problem  is  approached  using  Laplace  transform 
theory  for  the  time  behavior  and  Green  function  theory  for  the  spatial  structure.  Results  show 
that  the  growth  of  the  atmospheric  perturbations  closely  follows  the  growth  of  the  convectiv® 
boundary  layer;  the  transient  being  characterized  by  a  number  of  inertia-gravity  oscillations  of 
decreasing  intensity.  The  vertical  scale  is  determined  by  the  depth  of  the  convective  boundary 
layer:  and  the  horizontal  scale  is  determined  by  the  local  Rossby  deformation  radius.  Sinu¬ 
soidally  periodic  thermal  forcing  induce  periodic  atmospheric  cells  of  the  same  horizontal  scale. 
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The  intensity  of  mesoscale  cells  increases  for  increasing  values  of  the  wave  number,  reaches  its 
maximum  value  when  the  wavelength  of  the  forcing  is  of  the  order  of  the  local  Rossby  radius, 
and  then  decreases  as  the  wavelength  of  the  forcing  decreases.  However,  because  of  the  destruc¬ 
tive  interference  between  mesoscale  cells,  the  intensity  of  the  vertical  velocity  is  only  a  weak 
function  of  the  wave  numbers,  for  large  wave  numbers.  Periodic  square  wave  surface  thermal 
inhomogeneities  are  more  effective  than  sinusoidal  waves  in  generating  mesoscale  cells,  i.e.  the 
intensity  of  the  flow  is  generally  stronger. 


1  Introduction 

Subgrid-scale  parameterization  is  of  considerable  importance  in  regional  and  larger  scale  models 
because  these  models  seldom  have  more  than  few  grid  points  within  a  Rossby  radius,  yet  the 
mesoscale  flow  driven  by  thermal  horizontal  inhomogeneities  in  the  Convective  Boundary  Layer 
may  be  significant.  Therefore  mesoscale  effects  need  to  be  introduced  in  parametric  form  into  those 
large  scale  models. 

Terrain  in  the  real  world  is  almost  never  uniform,  for  instance,  in  natural  areas  there  may  be 
a  variety  of  vegetation  types  and  patches  of  vegetated  or  unvegetated  areas  which  can  modify  the 
mesoscale  atmospheric  flow.  Horizontal  thermal  inhomogeneities  in  the  planetary  boundary  layer 
are  a  well  known  source  of  mesoscale  circulation  systems.  Thermally-forced  mesoscale  systems  such 
as  land  and  sea  breezes,  mountain-valley  winds,  and  urban  heat  island  circulations  have  been  the 
subject  of  many  studies  (e.g.,  see  references  in  Pielke,  1984,  and  Pielke  and  Segal.  1986).  Some 
thermally-generated  mesoscale  circulations  are  due  to  horizontal  inhomogeneities  in  ground  wetness, 
vegetation  cover,  snow  cover,  cloud  cover  etc  (e.g.  Zhang  and  Anthes,  1982;  Garrett,  1982:  Van  and 
Anthes.  1988:  Segal,  et  al.,  1988).  In  addition  there  are  also  anthropogenic  modifications  of  soil 
and  vegetation,  such  as  irrigated  or  non-irrigated  fields,  crops  at  different  stages  of  development, 
inhabited  areas,  etc.  (Hanna  and  Sw'isher,  1971:  Anthes,  1984;  Segal  et  al.,  1984:  Segal  et  al..  1989). 
Such  thermal  inhomogeneities  may  extend  from  few  square  meters  to  hundreds  or  even  thousands 
of  square  kilometers  (Wetzel  and  Chang,  1988).  A  number  of  studies,  observational,  numerical  and 
analytical,  have  been  made  of  these  physiographically  forced  circulations  (Abbs  and  Pielke.  1986: 
Mahfouf  et  al.,  1987;  Mahrer  and  Pielke.  1977;  McCumber  and  Pielke,  1981;  Zhang  and  Anthes, 
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1982;  Yan  and  Anthes,  1988). 

Since  the  fine  horizontal  structure  of  the  terrain  thermal  inhoniogeneities  is  not  usually  resolved, 
however,  in  regional  and  large  scale  non-linear  numerical  models,  it  is  important  to  give  an  analyt¬ 
ical  quantitative  evaluation  of  their  impact  on  the  atmospheric  flow.  In  this  paper,  we  present  an 
analytical  theory  for  resolving  the  atmospheric  perturbation  and  for  evaluating  the  vertical  veloc¬ 
ities,  resulting  from  thermal  inhomogeneities  in  the  absence  of  large  scale  flow  (or  in  presence  of 
very  weak  large  scale  flow)  in  the  limit  of  a  linear  theory. 


2  Governing  Equation  and  Evolution  of  the  Horizontal  and  Ver¬ 
tical  Scales 


From  Rotunno  (1983),  and  Dalu  and  Pielke  (1989)  the  linearized  two-dimensional  primitive  equa¬ 
tions  describing  thermally-forced  atmospheric  flow,  in  the  absence  of  a  large-scale  wind,  can  be 
reduced  to  a  single  equation  for  the  streamfunction: 


where  /  is  the  Coriolis  parameter,  N  is  the  Brunt-Vaisala  frequency,  and  A"’  is  the  damping  time 
due  to  friction,  i.e.  the  persistency  of  the  mesoscale  flow  when  the  diabatic  forcing  vanishes.  In  this 
paper  we  use  the  following  values  for  the  Brunt-Vaisala  frequency,  inertia  frequency  and  friction 
time  scale: 


A’  =  10  ^  s  f  =  2u;sinQ  =  0(10  s  A  =  0(u;); 


12) 


Q  is  the  diabatic  buoyancy  source 


Q^_gfo]j^M!.He{h-z)  =  QoT{x)q{t)He(h-=) 


pcj, 


0 


(3) 


uniformly  distributed  through  the  depth  h  of  the  convective  layer  (CBL),  as  indicated  by  the 
Heaviside  function  He  in  eq.  (3).  FQT{x)q[i)\i  the  heal  flux  at  the  ground,  which  has  a  horizontal 
distribution  r(i)  and  a  time  behavior  q{t). 

(a)  Non-dimensional  Time  Variables  ,  .  ' 

'  '  j  iJl  I  lb 

fc  I - 


Statement  A  per  telecon 

Dr.  Robert  Abbey  ONR/Code  1122 

Arlington,  VA  22217-5000 

NWW  1/6/92 


VVe  first  define  the  characteristic  time  scale  of  the  system  T,  and  we  use  it  to  non-dimensionalize 
time: 


T  = 


=  O(10‘')s  r  =  tT-^ 


(4) 


T  decreases  with  increasing  latitude  and  friction  (friction  keeps  the  time  scale  T  finite  at  the 
equator).  We  introduce  the  Laplace  transform,  defined  as 

OC 

/(s)  =  I  {/(r))  =  y /(r)exp(-sr)dT  and  IX'’ {/(r)}  = /(r) 

0 

We  recall  also  that: 


/(r  — »  oc)  =  lim  sf{s) 

5— U 


(5) 


i.e.  /(r)  s=  sf{s)  when  s  <<  1  and  t  =  tq  >>  1  =>  t  =  to  >>  T  (6) 

to  is  the  time  when  an  almost  stationary  state  is  reached;  to  is  larger  than  the  characteristic  response 
time,  T,  of  the  system. 

(b)  Evolution  of  the  Vertical  and  Horizontal  Scales,  and  the  Non-dimensional  Coordinates 
We  recall  the  relation  between  the  diabatic  forcing  Q  and  the  height  h  of  the  CBL  (Green  and 
Dalu.  1980): 


r{x)  h{t)  =  rix)^^  dt' q{ty,  ho=^^J^  dt' q(t')  (7) 

The  Laplace  transforms  of  the  CBL  height,  h  and  the  non-dimensional  vertical  coordinate,  f].  are: 

d('S)  ' 

h  =  ho—;fi=-r  (8) 

3  h 

The  Laplace  transforms  of  the  horizontal  scale,  R,  and  the  non-dimensional  horizontal  coordinate, 
are: 


V/’  +  pV  <JP*p‘ 


p  =  3  -h  A; 


The  horizontal  scale  evolves  in  time  as: 

r‘  dt'  h(t  -  t') 


R{t) 


=  Ro  f 
Jo 


ho 


exp(-A/')  Joift'y,  Ro  =  ho 


N 

v/T^TF 


(9) 


(10) 
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Friction  acts  to  reduce  the  horizontal  scale  Rq  and  exponentially  damps  the  inertial  oscillations 
described  by  the  Bessel  function,  which  provides  some  damping  even  with  no  friction,  eq.  (10). 

The  evolution  of  the  vertical  scale  h(t)  and  of  the  horizontal  scale  R{t)  are  shown  in  Fig.la 
for  a  sinusoidal  in  time  thermal  forcing,  and  Fig. lb  for  an  impulsive  thermal  forcing;  for  different 
values  of  friction.  In  both  cases,  the  growing  time  of  the  boundary  layer  is  about  half  a  day, 
to  as  day.  The  horizontal  scale  grows  as  the  convective  boundary  develops,  through  inertia-gravity 
oscillations.  The  relative  amplitude  of  these  wave  is  of  the  order  of  0[exp(-At)  Joift)]:  i.e.  the 
inertia-gravity  oscillations  decay  as  a  Bessel  function  or  faster  (exponentially)  in  the  presence  of 
friction,  and  their  amplitude  is  usually  negligible  when  t  >  to  ^  0.5  day. 

(c)  Local  Rossby  Radius 

When  two  adjacent  regions  both  experience  a  warming  of  similar  intensity,  the  diabatic  heat,  in 
eq.(3),  is  given  by 

Q  =  QQq(t)[He{x)  +  (I  -  'r)He(-x)]He{ho-  z)  (11) 

with  7  <  0(1)  indicating  that  the  two  regions  have  similar  heating;  then  the  Brunt- V’aisala  fre¬ 
quency  difference  between  the  two  regions  is  of  the  order  of 

=  (12) 


The  Rossby  radius  Ro  has  been  replaced  by  the  local  Rossby  radius  Ro^  equal  to 

1/2 


=  ^ho 


(13) 


Because  of  the  reduction  of  the  Brunt-Vaisala  frequency  (eq.  12\  the  horizontal  scale  (eq.  13) 
is  correspondingly  reduced,  while  the  vertical  scale  remains  of  the  order  of  h(t).  The  aspect  ratio 
becomes: 


A  = 


I 


Atmospheric  cells,  generated  by  thermal  inhomogeneities  having  a  wavelength  smaller  than  twice 
the  local  Rossby  radius,  will  interact  with  each  other. 

(d)  Non-dimensional  Equation  for  the  Streamfunction 
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(14) 


Using  the  definitions  stated  above,  the  streamfunction  equation  (1)  can  be  written  as: 


-  d-  ■  -  d  - 

—ru  +  =  -Q(p)1Lq 


U'^  +  pV  \P+p^) 


\/p+p^ 


(15) 


The  variables  with  the  tilde  are  non-dimensional;  the  variables  with  the  hat  are  their  Laplace 
transforms. 

(e)  Boundary  Conditions  and  Green  Functions 

We  assume  that  the  terrain  is  flat  and  that  the  verticjil  velocity  vanishes  at  the  ground: 


w{x,z  -0,t)  =  w{i,fi  =  0,r)  =  w{i,f)  =  0,s)  =  0 


When  the  operator  in  eq.  (14)  is  elliptic,  the  Green  function  is: 


9<i>  =  -rr-in 


1/2 


2^^  L(^-C)"  +  (^-^02J 

The  perturbation  is  trapped,  i.e.  vanishes  at  oo: 


(16) 


(17) 


ti'(^  =  oo,q  =  oo,.«i)  =  0 


(18) 


When  the  operator  in  eq.  (14)  is  hyperbolic,  the  Green  function  is: 

g,.  =  i  [He{fi  -f,'-  \i  -  i'\)  -  Heifi  +  f}' -  -  ^'1)]  (19) 

The  perturbation  is  a  propagating  wave. 

Important  Note: 

When  the  operator  is  elliptic,  the  intensity  of  the  perturbation  decays  as  the  logarithm  of  the 
ratio  between  distance  from  the  source  and  from  its  mirror  image  (^',-fi'),  eq.  (17)  (the 

perturbation  is  confined  in  the  neighborhood  of  the  source  and  of  its  image).  The  presence  of  poles 
in  inverting  the  solution  in  the  transformed  space,  gives  propagating  inertia-gravity  waves  vanishing 
at  oo  in  the  physical  space;  the  related  solutions  are  in  terms  of  Bessel  and  trigonometric  functions 
(Dalu  and  Pielke,  1989). 


The  case  when  the  operator  becomes  hyperbolic  (non-vanishing  propagating  waves)  under  sinu¬ 
soidal  forcing  has  been  treated  by  Rotunno  (1983),  and  by  Dalu  and  Pielke  (1989).  Furthermore. 
Sun  and  Orlanski  (1981)  have  shown  that  some  modes  can  become  unstable  for  some  values  of  the 
pulsation  and  the  wave  number. 


3  Atmospheric  Response 

The  atmospheric  response  is  given  by  the  inverse  Laplace  transform  of  eq.  (20) 

V  =  Qoq(s)3ip)H,^ii,  fj)-,  ii  =  QQqis)3(p)Hu{t  q):  w  =  Qoqis)Mp)Hji,  q); 


f  - 

t  =  -  -u; 

P 


i  =  i  [g  -  A-^u 


(20) 


...  fh  .  . 

where  G(^  -  q)  =  /  dq'g(^  -('^q-  q') 
Jo 


The  integration  limit,  h,  is  the  Laplace  transform  of  the  depth  of  the  CBL.  In  addition. 


H{Lq)  =  J^dlG(i-i',q)^ 


D  is  the  horizontal  region  with  non-zero  forcing.  When  D  is  finite,  H  can  be  written  also  as: 


H(lq)  =  - 


Where 


Gu{^,q)  =  ^ 


1/2 

1/2 

e  +  q\ 

-  2  tan  ’  1 

V  W 

Ki)  ' 

<  i  )\ 

tan  1  — : — 

1  -  2  tan”’ 

-1-  tan”’ 

V  i  } 

\U 

V  ^  )\ 

+iv+  l)ln  + 


+(7)-l)ln[^n(^-l)']'^'|  (21) 

In  physical  units  we  have: 

0(x,2,O=  Y^{9(f)*  {exp(-At)-{yo(/0*-fft^(i.2.0}}}  (22) 

w'here  {  ♦  }  denotes  the  convolution  product  in  time.  The  time-dependent  transfer  functions 

can  be  written  in  terms  of  Bessel  functions,  trigonometric  functions  and  exponentials  (see  Dalu  and 
Pielke,  1989).  The  time  evolution  is  then  reconstructed  through  the  use  of  the  Faltung  theorem 
(Fodor,  1965).  The  resulting  behavior  is  similar  to  the  one  shown  in  Fig.l  for  the  Rossby  radius; 
i.e.  the  flow  evolves  through  a  series  of  damped  inertial-gravity  waves,  the  oscillations  are  not 
generally  in-phase  at  different  locations.  Making  use  of  the  knowledge  that  the  flow  decays  almost 
exponentially,  i.e.  with  an  e-folding  distance  equal  to  the  Rossby  radius,  as  shown  in  Fig. 2,  the 
structure  function  can  be  approximated  as: 


fj)  as  G,i,{i  =  0,f))  exp  (-Ifl) 

Stationary  State:  When  <>to  and  r(x)  =  H€(x),  we  have: 

I  /->  If  \  -,1.  /  Qo^o^o  Nhl 

V-' =  tl'oG,i,(^,T])  with  00= — r —  =  — 

A  2 


(23) 


(24) 


The  flow  from  eq.  (24),  shown  in  Fig.3a,  is  approximately  contained  in  a  region  2Ro  by  2ho.  The 
horizontal  and  the  vertical  momentum  components  are: 


u  =  ‘UoGm,t])  with  uo  =  —  =  —^=  — 

tlQ  Z  Zi 


^  ho  00  ho 

W  =  U'o  GJi.  Tj)  with  U:o  =  Uo—  =  —  =  — 


(25) 

(26) 


with  Tj  =  — ; 

ho 


^  Ro 


From  eq.  (25)  and  eq.  (26)  we  deduce  that: 
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i.e.  the  advection  terms  are  smaller  than  the  linear  term,  but  of  the  same  order;  thus  the  linear 
theory  underestimates  the  perturbations  predicted  by  a  nonlinear  theory.  For  a  quantitative  evalu¬ 
ation  of  the  nonlinear  contribution  see  the  paper  by  Pielke  et  a/ (1991).  The  along- front  momentum 
component  is: 

(27) 

The  temperature  perturbation  is: 

9  =  H€(x)He(ho-  z)  {ho  -  z)Q,  -  [Hei-x)  +  He(x)H€(z  -  he)]  (28) 

Here  we  refer  to  heated  land  but  the  temperature  perturbation  will  be  positive  whenever  a  CBL 
occurs  due  to  superadiabatic  lapse  rate  at  the  surface  (e.q.  over  warm  water  patches  also).  A:  is 
the  vertical  displacement  of  the  air  parcels: 

A^(T,^,to)  *  ~uix,z,to) 

When  the  flow  results  from  adjacent  differentially  heated  regions  in  eq.  (11),  the  intensity  of  the 
mesoscale  flow  reduces  to: 


=  7  Vb  (29) 

The  flow  intensity  is  proportional  to  7.  If  we  set  the  threshold  for  the  generation  of  significant 
mesoscale  flow  at  20  %,  the  influence  of  the  mesoscale  flow  need  to  be  included  (explicitely  or  in 
parametric  form)  when: 

^  >  0.20  =>  7  >  0.20 

Vo 

We  will  refer  to  this  measure  as  the  resolution  threshold  suggesting  that  for  7  <  0.2.  the  influence 
of  the  spatial  hetereogenity  can  be  ignored.  The  flow  shown  in  Fig. 3b  refers  to  the  case  when 
7  =  0.5  in  eq.  (11),  i.e.  the  upstream  region  experiences  50  %  less  heat  than  downstream  of  the 
thermal  discontinuity.  The  mesoscale  flow  shown  in  Fig.3b  is  less  intense,  (from  ,  eq.  (13)),  and  is 
less  extensive  than  the  flow  shown  in  Fig.3a. 
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4  Atmospheric  Response  to  Periodic  Thermal  Forcing 


Throughout  this  section  in  evaluating  the  averages  and  the  convolution  integrals  along  the  horizontal 
coordinate  we  make  use  of  the  fact  that  the  intensity  of  the  flow  decays  exponentially  with  an  e- 
folding  distance  equal  to  the  Rossby  radius  (Fig.2); 


g{z)  =  G^.{x 


2;r  I  z  -  /in 


When  the  diabatic  heat  is  sinusoidally  periodic  over  the  entire  horizontal  domain,  with  as 
wavelength: 

Q  =  Qoq(t)  He(h  +  \  sin  ;  m  =  ^  (31) 

then  the  streamfunction  is  given  by 

,  /  N  fm‘Kx\ 

Here  m  is  the  number  of  wavelengths  in  a  2Rq  distance  (m  need  not  be  integer).  If  we  set  the 


resolution  threshold  at  0.20: 

V’o 

^  .w  n  rtf 


>  0.20  =>  - - rr — -  >  0.20 

(mjr)2  +  1 


0.4  <  m  <  1.2 


Thus  the  mesoscale  flow  need  to  be  included  (explicitely  or  in  parametric  form)  when  0.4  <  m  <  1.2. 
The  intensity  of  the  atmospheric  response  as  a  function  of  m  is  shown  in  Fig.4a.  The  intensity 
of  the  mesoscale  cells  monotonically  decreases  for  m  large  (see  also  shown  in  Fig. 5).  The  vertical 
velocity  is  : 

,  (mff)2  .  /m7rx\ 


dtl’  ,  ,  (rriTry 


.  /  Tmrx\ 
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When  m  >  1,  the  vertical  velocity  is  a  weak  function  of  the  wave  number  (Fig.4c).  When  the 
forcing  is  a  periodic  square  wave  over  the  entire  horizontal  domain,  with  Lm  as  wavelength: 


Q  =  Qoq{t)  He{h  -  z) 
then  the  streamfunction  is  given  by 


1  1  ^  4  (mn-Kx\ 

5  +  5  ^ 


n=1.3-- 


V  Ro  ) 


m  = 


2Ro 


(34) 


v=il'og{z)  ^ 


mnn 


r,^,3  T2T(mnff)2  +  1 

If  we  set  the  resolution  threshold  at  0.20; 

4  mnir 


cos 


(  mn'Kx\ 

V  Rq  } 


(35) 


>  0.20 


>  0.20 


0.05  <  m  <  2.2 


Thus  for  square  heating  patches,  the  mesoscale  flow  needs  to  be  included  (explicitely  or  in  para¬ 
metric  form)  when  0.05  <  m  <  2.2.  The  intensity  of  flow  as  a  function  of  m  for  a  periodic  square 
forcing  is  shown  in  Fig. 4b,  and  the  circulation  cells  are  shown  in  Fig.6.  The  intensity  of  the  flow, 
in  general,  is  larger  for  a  periodic  square  forcing  than  for  a  sinusoidal  forcing;  the  difference  being 
larger  at  small  wave  number.  At  large  wave  number  the  mesoscale  cells  interact  destructively, 
reducing  the  intensity  of  the  flow.  The  vertical  velocity  (Fig.4d)  is  : 


w 


di'  /  X  4  (mnTr)^  .  f mm:x\ 

=  -—  =  wog(2)  >  — - — r — 

dx  nx  (mn7r)2 -f  1  \  Ro  J 


n=1.3- 


(36) 


and  the  temperature  perturbation  is: 

e  =  Heiho  -  :) 


1 

1 

OC 

4  .  - 

f  mn'Kx\ 

2 

+  2  , 

E 

n  =  1.3-- 

—  sin 
m:  ' 

K  Ro  ) 

(ho  -  z)Q, 


I 


-o- 


He(ho  -  z) 


1  1  ^  4  /  mmrx\  1 

5  +  JJI 


wQ,j 


(37) 


5  Conclusions 


Making  use  of  the  theory  developed  by  Dalu  and  Pielke  (1989),  we  study  the  atmospheric  response 
to  spatial  gradients  of  diabatic  heat  fluxes  in  the  convective  boundary  layer.  Results  show  that, 
for  sudden  changes,  the  atmospheric  perturbation  decays  exponentially  with  an  « -folding  distance 
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equal  to  the  local  Rossby  radius.  Square  wave  periodic  surface  thermal  inhomogeneities  generate 
stronger  atmospheric  cells  than  sinusoidally  distributed  inhomogeneilies.  The  intensity  of  the  flow 
increases  as  the  wave  number  increases,  reaches  its  maximum  when  the  wave  length  is  of  the  order 
of  the  Rossby  radius,  and  then  decreases  linearly  for  large  wave  number.  Thus  significant  atmo¬ 
spheric  cells  are  generated  by  periodic  thermal  inhomogeneous  patches  only  when  the  wavelength 
is  a  significant  fraction  of  the  local  Rossby  radius.  It  la  unnecessary  therefore  to  explicitly  resolve 
in  num  cal  models  thermal  inhomogeneities  on  scales  much  smaller  than  the  local  Rossby  radius 
(e.g.  as  assumed  in  Avissar  and  Pielke,  1989).  However,  since  at  high  wave  numbers,  the  ver¬ 
tical  velocities  is  significant  and  almost  constant,  mesoscale  heat  fluxes  may  be  significant;  these 
subgrid-scale  effects  need  to  be  introduced  in  parametric  form  in  LAM  and  GCM,  which  have 
only  few  grid  points  within  a  Rossby  radius.  Solutions  for  the  atmospheric  perturbations  in  the 
presence  of  light  ambient  winds,  can  be  derived  from  the  theory  presented  in  this  papei,  and  an 
extension  of  our  theory  in  this  direction  is  desirable.  When  the  ambient  flow  is  strong,  the  pattern 
of  the  perturbations  is  quite  different  (Hsu,  1987;  Robichaud  and  Lin,  1989)  however  the  intensity 
of  the  meso'  cale  perturbation  s  negligible,  since  it  decays  as  l/U^. 
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Figure  Captions 


Fig.l  The  evolution  of  the  Rossby  radius  for  different  values  of  friction:  -a-  for  a  sinusoidal 
forcing;  h{t)  =  ho(l  —  cosiJt)/2  is  given  in  the  insert;  -b-  for  a  convective  boundary  layer  asymp¬ 
totically  growing  to  /?o;  h(t)  =  /io(l  -  exp(-//<o)  is  given  in  the  insert. 

Fig. 2  -a-  Streamfunction.  v  [10^  m^/s),  at  z  =  ho  versus  the  horizontal  distance  from  the  thermal 
discontinuity.  The  function  v(J  =  0.0, r  =  ho)  exp(— |i|/ro)  is  plotted  to  show  that  the  intensity 
of  the  flow  decays  essentially  e.xponentially  with  an  e-folding  distance  equal  to  the  Rossby  radius 
(broken  line),  -b-  Vertical  profile  of  v'  at  z-  =  0-0. 

Fig. 3  -a-  Streamfunction  isolines  for  a  discontinuous  change  of  heat  flux  at  i  =  0.0  at  40°  of 
latitude  (A  =  u;);  hy  =  1000  m.  At’  =  0.2  [10^  m^/s].  -b-  Streamfunction  isolines  w'hen  one  region 
experiences  a  heat  flux  50  less  intense  than  the  adjacent  region,  i.e.  y  =  0.5. 

Fig. 4  -a-  Intensity  of  the  flow,  induced  by  a  sinusoidal  periodic  forcing,  versus  horizontal  w’ave  num¬ 
ber;  -b-  same  as  -a-  but  for  a  periodic  square  wave;  -c-  vertical  velocity  for  a  sinusoidal  periodic 
forcing;  -d-  vertical  velocity  for  a  periodic  square  wave;  -e-  average  vertical  mesoscale  heat  flux 
for  a  sinusoidal  periodic  forcing;  -f-  average  vertical  mesoscale  heat  flux  for  a  periodic  square  wave. 

Fig. 5.  -  Mesoscale  atmospheric  cells  for  four  values  of  the  horizontal  wave  number  for  a  sinu¬ 
soidal  periodic  forcing.  Af  =  0.2  [10^  m^/s]. 

Fig. 6.  -  Mesoscale  atmospheric  cells  for  four  values  of  the  horizontal  wave  number  for  a  square 
wave  periodic  forcing,  Ac  =  0.2  [10^  m^/s]. 
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